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Abstract 

We present a complete computation of the inflaton decay into very massive 
fermions during preheating. We show that heavy fermions are produced very 
efficiently up to masses of order 10^^-10^^ GeV; the accessible mass range 
is thus even broader than the one for heavy bosons. We apply our find- 
ings to the leptogenesis scenario, proposing a new version of it, in which the 
massive right-handed neutrinos, responsible for the generation of the baryon 
asymmetry, are produced during preheating. We also discuss other produc- 
tion mechanisms of right-handed neutrinos in the early Universe, identifying 
the neutrino mass parameters for which the observed baryon asymmetry is 
reproduced. 
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1 Introduction 



It is commonly believed that the Universe underwent an early era of cosmological infla- 
tion |]T| . The flatness and the horizon problems of the standard big bang cosmology are 
indeed elegantly solved if, during the evolution of the early Universe, the energy density 
happens to be dominated by the vacuum energy of a scalar field - the inflaton - and 
comoving scales grow quasi-exponentially. 

At the end of inflation the Universe was in a cold, low-entropy state with few degrees 
of freedom, very much unlike the present hot, high-entropy Universe. At this stage, the 
Universe does not contain any matter and therefore it looks perfectly baryon symmetric. 
However, considerations about how the light element abundances were formed when the 
Universe was about 1 MeV hot lead us to conclude that ub/s = (2-9) x 10^^^. Here 
ub/ s is the difference between the number density of baryons and that of antibaryons, 
normalized to the entropy density of the Universe. 

Until now, several mechanisms for the generation of the baryon {B) asymmetry 
have been proposed 0. Grand Unified Theories (GUTs) unify the strong with the 
electroweak forces and predict baryon-number violating interactions at tree level. In 
these theories, the out-of-equilibrium decay of heavy Higgs particles can indeed explain 
the observed baryon asymmetry. In the theory of electroweak baryogenesis, baryon 
number violation takes place at the quantum level, caused by unsuppressed baryon- 
number violating sphaleron transitions in the hot plasma ||^. 

Since B and L - where L is the lepton number - are reprocessed by sphaleron 
transitions, while the anomaly-free linear combination B — L is left unchanged, the 
baryon asymmetry may be generated from a lepton asymmetry Indeed, once the 
lepton number is produced, thermal scatterings redistribute the charges and convert 
(a fraction of) L into baryon number. In the high-temperature phase of the standard 
model, the asymmetries of baryon number B and of B — L are therefore proportional ||^: 

22ng + 13nH) ' 

where riH is the number of Higgs doublets and Ug is the number of fermion generations. 

In the standard model as well as in its unified extension based on the group SU{5), 
B — L is conserved and no asymmetry in B — L can be generated. However, adding 
massive right-handed Majorana neutrinos to the standard model breaks B — L and 



B = a{B-L), 
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the primordial lepton asymmetry may be generated by their out-of-equihbrium decay. 
This simple extension of the standard model can be embedded into GUTs, as in the 
case of 50(10). Heavy right-handed Majorana neutrinos are the key ingredient to 
explain the smallness of the light neutrino masses via the see-saw mechanism 0. The 
presence of neutrino masses and mixings seems to be the most natural explanation of the 
recent reports from the Super-Kamiokande [|7| and other collaborations indicating 
the existence of neutrino oscillations. In light of these considerations, the generation of 
the baryon asymmetry through leptogenesis looks particularly attractive. 

The leptogenesis scenario depends crucially on the mechanism that was responsible 
for populating the early Universe with right-handed neutrinos, and consequently on the 
thermal history of the Universe, and on the fine details of the reheating process after 
inflation. One goal of this paper is to discuss several production mechanisms of heavy 
right-handed neutrinos in the Universe, to compare them and to identify the regions of 
the appropriate parameter space where the production mechanism is efficient enough 
to explain the observed baryon asymmetry. 

The simplest way to envision the reheating process after inflation is to assume that 
the comoving energy density in the zero mode of the infiaton decays perturbatively into 
ordinary particles, which then scatter to form a thermal background 0. It is usually 
assumed that the decay width of this process is the same as the decay width of a free 
inflaton field. Of particular interest is a quantity known as the reheat temperature, 
denoted as Tun- This is calculated by assuming an instantaneous conversion of the 
energy density in the infiaton field into radiation when the decay width of the infiaton 
r,^ is equal to H, the expansion rate of the Universe. This yields 



■RH 



r^Mp, (2) 



where Mp is the Planck mass. The commonly-accepted assumption is that the heavy 
right-handed neutrinos with mass Mn were as abundant as photons at very high tem- 
peratures. This assumption requires not only that Trh - Mn, but also that the heavy 
neutrinos are abundantly produced by thermal scatterings during the reheating stage. 
This condition, as we will discuss in sect. 3.1, significantly limits the allowed range of 
neutrino masses compatible with leptogenesis. 

There might be one more problem associated with the hypothesis that T^h ^ Mn 
in the old theory of reheating, and that is the problem of relic gravitinos [110]. If one 
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has to invoke supersymmetry to preserve the flatness of the inflaton potential, it is 
mandatory to consider the cosmological implications of the gravitino - the spin-3/2 
partner of the graviton which appears in the extension of global supersymmetry to 
supergravity. The slow gravitino decay rate leads to a cosmological problem because the 
decay products of the gravitino destroy light nuclei by photodissociation and hadronic 
showers, thus ruining the successful predictions of nucleosynthesis. The requirement 
that not too many gravitinos are produced after inflation provides an upper bound on 
the reheating temperature Tjih of about 10*^-10^° GeV, depending on the value of the 
gravitino mass |[Tl|]. In the following, T^h will be therefore intended as the largest 



temperature allowed after inflation from considerations of the gravitino problem. 

In order to relax the limit on Mn imposed by the gravitino problem, we will consider 
the possibility that the heavy neutrinos are produced directly through the inflaton decay 
process. This is kinematically accessible whenever 

Mn < m^, (3) 

where is the inflaton mass. In the case of chaotic inflation with quadratic potential, 
the density and temperature fluctuations observed in the present Universe determine 
and require to be smaller than about 10^^ GeV. 

The outlook for leptogenesis might be brightened even further with the realization 



that reheating may differ significantly from the simple picture described above [|12|, 



1^ , [15| , p!q] . In the first stage of reheating, called preheating [0, nonlinear quantum 
effects may lead to extremely effective dissipative dynamics and explosive particle pro- 
duction, even when single particle decay is kinematically forbidden. In this picture, 
particles can be produced in a regime of broad parametric resonance, and it is possible 
that a significant fraction of the energy stored in the form of coherent inflaton oscilla- 
tions at the end of inflation is released after only a dozen or so oscillation periods of the 
inflaton. The preheating stage occurs because, for some parameter ranges, there are new 
non-perturbative decay channels. In the case of bosonic particle production, coherent 
oscillations of the inflaton field induce stimulated particle emissions into energy bands 
with large occupancy numbers [|I2|. The modes in these bands can be understood as 
Bose condensates, and they behave like classical waves. A crucial observation for GUT 
baryogenesis induced by the decay of very heavy Higgs bosons is that particles with 
masses larger than that of the inflaton may be easily produced during preheating . 
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Indeed, for coupling constants of order unity one would have copious production of bo- 
son particles as heavy as 10^^ GeV, i.e. 100 times greater than the inflaton mass. This 
is a major departure from the old constraint of reheating. 

In this paper we wish to provide the first complete calculation of the infiaton decay 
into heavy fermions during preheating. We will show that fermion production is ex- 
tremely efficient, in a mass range even broader than the one for heavy bosons: fermions 
may be generated up to masses of order of 10^^ GeV. What distinguishes the production 
of very massive fermions and bosons in an oscillating background is the expression for 
the total mass. For bosons, the total mass can never vanish and the production reaches 
the maximum when the amplitude of the inflaton goes through zero. For fermions, 
the total mass can vanish for particular values of the inflaton fleld, rendering particle 
creation much easier. 

We will numerically compute the density of the massive fermions produced at the res- 
onance stage. This is the crucial parameter for Icptogenesis, when these heavy fermions 
are identified with the right-handed neutrinos giving rise to the lepton asymmetry. We 
want to stress that the out-of-equilibrium condition is naturally achieved in this sce- 
nario, since the distribution function of the fermionic quanta generated at the resonance 
is far from a thermal distribution. We will show that the observed baryon asymmetry 
may be explained by the phenomenon of leptogenesis after preheating, with a reheating 
temperature compatible with the gravitino problem. 

The paper is organized as follows. In sect. 2 we present our calculation, with 
numerical results as well as some analytical estimates concerning the production of 
massive fermions during preheating. In sect. 3 we discuss and compare the relevant 
production mechanisms of right-handed neutrinos in the early Universe, identifying the 
appropriate neutrino-mass parameters where these mechanisms are the most efficient, 
as far as the generation of the baryon asymmetry is concerned. 

2 Heavy Fermion Production at Preheating 

In this section we describe the basic physics underlying the mechanism of heavy fermion 
generation during the preheating stage, perform the relevant numerical calculations and 
present some analytical estimates. In the following we will focus on the model of chaotic 
inflation, with a massive inflaton having quadratic potential = |m^0^. Here 
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~ 10^^ GeV is fixed by the COBE normalization of tlie cosmic microwave background 
anisotropy. 

We will suppose that the infiaton field is coupled to a very massive Dirac fermion 
X with bare mass mx via the Yukawa coupling 

Cy = 9(j)XX. (4) 

The total mass of the fermion X is then given by 

m{t) = mx +g(t){t)- (5) 

When we apply our results to the leptogenesis scenario, the fermion X will be identified 
with the lightest of the right-handed Majorana neutrinos A^^i. Although in this section 
we will restrict our considerations to the case of Dirac particles, the treatment of Ma- 
jorana fermions is completely analogous. Our final results regarding the abundances of 
particles (with equal amount of anti-particles being produced) are valid for both Dirac 
and Majorana fermions. 

2.1 The Basic Formalism 

We start (see e.g. discussion in ref. [|l8l) by canonically quantizing the action of the 



massive field X in curved space with Friedmann-Robertson- Walker metric. In the sys- 
tem of coordinates in which the line element is given by ds^ = a^{ri){drf — dx'^), where 
a is the scale factor of the expanding Universe and rj is the conformal time defined as 
drj = dt/a, the Dirac equation becomes 

Here H = [a' /a?) is the Hubble rate, the prime denotes derivative with respect to 
conformal time, and the 7-matrices are defined in flat space-time. By defining x = 
a~^/^X, eq. (^) can be reduced to the more familiar form 

(i7'^9^-am)x = 0. (7) 

Since a is a function of 77, but not of x, spatial translations are symmetries of space- 
time, and we can separate the variables using the decomposition 

X{x) = J j^^e-'''-'^J2 [urik,v)arik) + Vrik,7])bli-k)] , (8) 
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where the summation is over spin, and Vr{k) = CuJ{—k). We impose the canonical 
anticommutation relations on the creation and annihilation operators 



{ar{k),al{k')} = {br{k),bl{k')} = Srs6{k-k') 



(9) 



which, together with the quantization conditions, determine the normalization of the 
spinors u, 



ul{k, r])usik, 7]) = vl{k, T])vs{k, t]) = 5rs, ul{k, T])v,{k, t]) = 0. 



(10) 



Equations (|iy) are valid at any conformal time, since they are preserved by the evolution. 

In the representation in which 7° ~ g ^^"^ with the definition u = ^ ' 

the equation of motion can be written as a set of uncoupled second-order differential 
equations 

r W2 _ 1 

u±{k) = (11 



dr] 



2 + cij^ ± i{a'm + am') 



LU 



We can now write the Hamiltonian as 



k-^ + m^a\ 



H{r])=J rfW Hdo) X = J d'kY,{Ek (r/) [aj(fc)a,(fc) - br{k)bl{k) 

+ Fk iv) bri-k)ar{k) + Fl (r]) al{k)bl{~k)} 

By using the equations of motion, we find 



(12) 



(13) 



Ek = 2k Re(?i+u_) + am (l - 2 
Fk = k (u"^ — + 2am u+u 



U.U4 



El + \Fk\^ = u:\ 



(14) 
(15) 
(16) 



Here we have chosen the momentum k along the third axis, and selected the gamma- 
matrix representation in which 7^ ~ g ^ • 

In order to give a "quasi-particle" interpretation, we diagonalize the Hamiltonian in 
eq. ([I3D with a time-dependent Bogolyubov canonical transformation, and define the 
new creation and annihilation operators 



a{k,r]) = a{k,ri)a{k) + (3{k,7])b\—k) 
b^{k, r]) = -(3*{k, 7])a{k) + a*{k, r])b^{-k). 



(17) 
(18) 



6 



Imposing canonical anticommutation relations on the operators a and b, we find + 
= 1. For 

a _ Ek + uj 1^,2 _ \Fk\'^ _ uj-Ek 
P ' l'^' 2Lu{LU + Ek) 2uj ' ^ ' 

the normal-ordered Hamiltonian in terms of the "quasi-particle" operators is diagonal, 

H{r]) = [ d^kJ2^{v)[al{k)ar{k) + bl{k)br{k)]. (20) 

Next, we define a "quasi-particle" vacuum, such that a|0^) = &|0^) = 0. The total 
number of produced particles up to time rj (equal to the number of produced antipar- 
ticles) is given by the vacuum expectation value of the particle number operator A^, 



n{^) = (0,|iV|0,) = ^— - / dk e (21) 

The density of produced particles is then computed by integrating the equations of 
motion ( |Tl| ) with an initial condition at time rj = given byQ 



«±(0) = (iT^), u'^{0)=tku^{0)T^amu±{0). (22) 

This boundary condition corresponds to = u;, = at r/ = or, in other words, to 
an initial vanishing particle density. 

2.2 Numerical Results 



The equations of motion ([ll[) describe oscillators with time-varying complex frequency. 
If m is constant, the time dependence enters only through the scale factor a. The 



corresponding gravitational creation of heavy fermions was studied in ref. ||T9|]. However, 
of particular interest is the case in which m = m{t) is a periodic function of time. This 
is realized when the scalar field 0, coupled to X as in eq. @) is homogeneous and 
oscillates in time with frequency V^"(0). It is useful to write the equations of motion 
in terms of dimensionless variables. We introduce a dimensionless time r = m^rj, as 
well as a dimensionless field if = 0/0(0), so that the scalar field is normalized by the 
condition v^(0) = 1. We define 0(0) as the value of the infiaton field at the moment 



when its oscillations begin, 0(0) ~ 0.28Mp [|T^. With these redefinitions, the equation 



of motion for the background field ip does not contain any parameter, while the fermion 

^In practice, one has to solve only the evolution equation for since ii_ is determined by the 
equation m_ — (amu-\- — iu'_^)/k. 
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mass is measured in units of and the strength of the fermion couphng to the external 
background is determined by the dimensionless combination g(p{Q)/m^. For the sake of 
correspondence with the bosonic case, we introduce the parameter 

q = g^ct>\Q)/Aml. (23) 

The production of fermions by an external oscillating background in Minkowski 
space-time was studied in refs. pP] , pT| , The production of massless fermions in a 
inflaton model, studied in ref. p2|, can be reformulated into particle production in 



an expanding Universe by a simple conformal transformation. Moreover, only the case 
of moderate q {q < 100) has been previously considered. The parametric production of 
massive fermions in an expanding Universe has never been analyzed before. By analogy 
with the bosonic case |jT2|, [1^, |1^ , we expect that an efficient production of very massive 
fermions will require a very large value of q. 

The result of our numerical integration is best summarized in fig. ^ For different 
values of the fermion mass mx and of the q parameter, fig. shows px / p, the fraction 
of the infiaton energy density p which ends up in the fermionic energy density px- In 
our simulation, we have neglected the back reaction of the produced fermions in the 
evolution of the infiaton field. 

Fermion production is efficient up to a time at which ratio px/p freezes out. For 
fixed q, the larger mx is, the earlier this freeze-out occurs. In particular, near the cut-off 
of px/p at large mx, the production ceases just after a few oscillations of the inflaton 
fleld. In flg. 1^ we have plotted the flnal phase-space density of produced fermions in 
comoving volume for mx = 100 m^, q = 10^ and q = 10^. For q = 10^ this distribution 
is reached after twenty inflaton oscillations, while for q = 10^ it is reached just after the 
flrst inflaton oscillation. In flg. |l|, the curves for q = 10^ and q = 10^ correspond to an 
evolution up to 20 inflaton oscillations. For q = 10^ and mx < 100 m^, the freeze-out 
has not been reached in 20 oscillations and px/p would have grown further, had we 
integrated for longer times. This explains the slope of the curves in flg. [l| at small mx 
for q > 10^. In general, if we integrate up to the freeze-out, we flnd that the ratio px/p 
is almost constant with mx up to a cut-off value (mx)max much larger than m^. In this 
regime, the number density of fermions is inversely proportional to mx- 

The time evolution of the phase-space density of produced fermions in the cases 
in which more than one inflaton oscillation is required to reach the flnal distribution 
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Figure 1: The fraction of the energy density of produced fermions with respect to the 
total energy density, as a function of the fermion mass mx in units of the infiaton mass, 
for various values of q. 



is shown in figs. ^ and ^. For the parameters shown in fig. ^, the final distribution 
was reached after three inflaton oscillations, while in the case of fig. |^ the X-particle 
production continues up to the twentieth oscillation. 

Particle production at large q and large mx is similar to the "instant preheating" 
phenomenon discussed in ref. p3|. However, contrary to what is discussed in ref. p3[ 



where the production of heavy fermions was possible only through the coupling of the 
fermion field to other bosonic degrees of freedom produced by preheating, in our case 
the instant preheating takes place by coupling the fermion directly to the inflaton field, 
and therefore this realization is much less model-dependent. 

Let us summarize our numerical results in a form suitable for the analytical estimates 
we will perform in sect. 2.3. Because of the Fermi-statistics, n{k) = < 1. The 
maximum value of n{k) ~ 1 is reached rapidly at some k = fcmax- With time, if particle 
creation is still efficient, /^max grows at each oscillation. We can estimate the number 
density of the produced fermions as n = / d^kn{k) oc k'^^^. From our numerical results 
we observe that, at fixed mx, the freeze-out value of fcmax scales as /cmax oc g", with 
a slightly larger than 1/3; see figs. |l| and |^. On the other hand, fcmax has to scale 
as mx to explain why px/p is nearly independent of mx, at constant q. We also 
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Figure 2: The final phase-space density of produced particles for two values of the 
parameter q. The X-fermions are taken to be 100 times heavier than the inflaton. At 
q — 10^ the freeze-out of the particle production was reached after the first infiaton 
oscillation, while for g = 10^ it required twenty oscillations. 




Figure 3: The phase-space density of produced particles after n = 1 and n = 3 infiaton 
oscillations for q = 10® and for X-fermions 100 times heavier than the infiaton. The 
distribution at n — 3 coincides with the final distribution. 
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Figure 4: The phase-space density of produced particles after n = 1 and n = 20 infiaton 
oscillations for q = 10^ and for X-fermions 100 times heavier than the infiaton. The 
distribution at = 20 coincides with the final distribution. 

find numerically that the cut-off value (mx)max, i-e. the maximum mass of produced 
fermions, is about (mx)max — 

2.3 Analytical Estimates 

The goal of this section is to provide some analytical estimates that can help us to 
understand the numerical results presented above. 

Compared with the bosonic case, the novel feature of the fermionic production is 
that much heavier particles can be generated rather efficiently during the oscillations of 
the infiaton field, see fig. The key difference between the Bose and Fermi cases resides 
in the expression for the effective particle mass in an external oscillating background. 
In the Fermi case it is given by eq. (|^), m{t) = mx + g4'{t), while in the Bose case 
we have ni?{t) = m\ + 9'^4''^{t)- In both cases, the production of particles with masses 
larger than the infiaton mass {i.e. the frequency of the infiaton oscillations) requires 
large q. This, in turn, means that m{t) is typically large, and creation of X-particles is 
impossible at all times, except for very short periods when m(t) fluctuates around its 
minimum value. In the bosonic case m(t) can never be smaller than mx. However, in 
the fermionic case, m{t) can vanish as long as the amplitude of the infiaton field is large 
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enough, |0| > mxig- 

In the expanding Universe the amphtude 0o of the oscillating inflaton field decreases 
with time. It is reasonable to assume that the sharp cut-off in the particle production 
at large mx that we observe in fig. [I| corresponds to a situation in which the condition 
m(t) = cannot be satified even during the first oscillation of the infiaton field. To 
verify this assumption, let us write m(t), with the help of eq. (pSj), as 



(j) ^JquiA) 
m = mx + 2^m^^^ ~ mx + cos(27rA^) . (24) 

Here we have taken into account that cj) oc when the energy density of the Universe 
is dominated by the oscillating infiaton field, and we have denoted with N the number 
of oscillations of the infiaton field = m^t/{2T[). 

At large mx, the minimum of m is reached around = 1/2. Therefore, if 

2 

mx > m0- y/q, (25) 

71 

the total mass m never vanishe^. The value of (mx)max given in eq. (^) is already 
in good agreement with our numerical results. However, eq. (Bl) cannot be trusted at 
small A^. Actually, our numerical integration shows that the minimum of the infiaton 
amplitude 0o during the first oscillation is given by 



0.25. (26) 



0(0) 

This means that the cut-off value of the mass is at 

("^x)max ^ m^^, (27) 

in perfect agreement with the results presented in fig. ||. 

If mx < (mx)max5 m{t) can vanish more than once during the infiaton oscillations. 
In particular, in the case in which mx ^ m^, one can suppose that the production of 
particles continues until the amplitude of the oscillations drops below the critical value 

(0o)cnt = — = /^0(O). (28) 
g 2m 

At later times, when 0o < (0o)crit, the ratio px/ p '^^ frozen, since particle production has 

stopped. Let us now estimate analytically the final value of px / P- During the infiaton 

^We are considering here the case in which g(t>{Q) and mx have the same sign. If these two terms 
have a relative minus sign, particle production can be extended to even larger values of mx- 
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oscillations, the Fermi distribution function is rapidly saturated up to some maximum 
value of the momentum k, i.e ~ 1 for k < /cmax and it is zero otherwise. The value 

of fcmax increases at each inflaton oscillation until particle production stops. Therefore, 
what is relevant for the determination of the final value of px/p is ^max at the freeze-out, 
i.e. when eq. ( PBD is satisfied and particles are no longer generated. At this moment, 

P.-^, (29) 

where pmax is the maximum physical momentum at freeze-out, pmax = ^max/ct- 

We estimate Pmax as the maximum value of the physical momentum above which 
the evolution of the vacuum mode function is adiabatic. We choose as a trial mode 
function the following expression 



U^V) = fl - ^) X exp (-Z roodv'], (30) 



,2V Lj 

where nies = ma. Notice that this solution corresponds to a vacuum, since it gives 
P{k) = 0. This trial function solves to a good approximation the equation of motion if 
the following condition holds 



1 , rries 



2 V UJ J 



UJ 2 \ UJ \ UJ 



<^uj\ (31) 



where we have used the property uj' = nicsm'^s/uj. As we have explained above, 
for large q, particle production takes place during short intervals when m^g ~ 0, or 
cos(m0t) ~ —mx/g<pQi neglecting the expansion of the Universe during these short 
periods of particle creation. Under these circumstances, the condition ( ^I]) reduces to 

2pmx + ^ -m\\ml<^Ap\ (32) 

where we have used ~ A; = pa. As a result, at the end of the particle production 
epoch, see eq. (p8D, we get 

pLx - i^mxml. (33) 



To compare with the numerical results of sect. 2.2, it is convenient to express eq. (|3 
in terms of the comoving momentum, 

mx 



krn.. = ) . (34) 
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This reproduces the approximate scahng law fcmax oc q^^^m^^^ observed in the numerical 
results of sect. 2.2. Notice that, at earlier epochs when (yf^0Q ^ m\, the value of Pmax 
is larger, but the corresponding fraction of X-particles produced at this stage becomes 
subdominant (because it is red-shifted) with respect to the fraction of particles generated 
right before the freeze-out. In terms of the parameter gefT(^) = 5'^0^(^)/4^«|, we obtain 
that at early times pmax oc g^g^, while at the end of the particle production epoch 

1/6 

J'max . 

Our result for Pmax amd (mx)max and their relative scaling with q for the fermions 
produced during preheating are different from what argued in ref. and from the 
behaviour of the same quantities encountered when dealing with boson particle pro- 
duction during preheating. In particular, for bosons, both Pmax and (mjs:)max scale like 

iL' m m- 



We can now use eqs. (^) and (0) to find that px — m|^m^/67r^. The inflaton 
energy density at freeze-out is p = 7ii?^(fP'/2 = 0^(O)m^/8g, and therefore the energy 
density fraction is 

^ ~ —^a -—a' (35) 
p -37r2 02(o)^-3vr2^- ^^^^ 

This expression describes quite well the behaviour observed in fig. |^ for a large range of 
q. Equation (|35|) does not depend upon rrix, it is approximately proportional to q and 



it gives a reasonable estimate for the overall magnitude. Indeed, for 0(0) ^ 0.28Mp [|14 



and m^/Mp ~ 10~^, eq. ( pSj) reduces to px/p — 2 x 10"^^g, in good agreement with 
our numerical results. At very large values of q the expression ( p5| ) underestimates the 
ratio px/p] for instance at g = 10^, it gives a value of px/p smaller than the numerical 
result by about one order of magnitude. 

In the next section, we will apply our results on fermionic preheating to the case of 
heavy right-handed neutrinos and leptogenesis. We want to emphasize, however, that 
the results obtained in this section may be relevant in other cosmological contexts, such 



as the generation of superheavy dark matter after inflation [£4 
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3 Leptogenesis 



The Lagrangian terms relevant for leptogenesis describe the interactions between the 
massive right-handed neutrinos A^, the lepton doublet ii, and the Higgs doublet H, 



C = -NKHiL - -N^MN + h.c. 



(36) 



Here the Yukawa couplings hi, and the Majorana mass M are 3x3 matrices and gener- 
ation indices are understood. We choose to work in a field basis in which M is diagonal 
with real and positive eigenvalues ordered increasingly (Mi < M2 < M3). The mass 
matrix of the light, nearly left-handed, neutrinos is given by 



= -hlM~^K{H)\ 



(37) 



The decays of the heavy neutrinos into leptons and Higgs bosons violate lepton 
number 



N 
N 



Hi, 
Hi. 



{3t 



The interference between the tree-level decay amplitude and the absorptive part of 
the one-loop diagram can lead to a lepton asymmetry of the right order of magnitude 
to explain the observed baryon asymmetry, as has been extensively discussed in the 
literature 0, |2^, ^ (for reviews, see ref. [^). The interference with the one-loop 
vertex amplitude yields a CP-violating decay asymmetry for A^i equal to 



1 



5:1m 

j=2,3 



Khl 



fix) 



1 - (1 x) In 



l + x 

X 



(39) 



(40) 



The absorptive part of the one-loop self-energy gives a contribution to the A'^i asymmetry 
which, in the case of only two-generation mixing, is given by 



Im 






2 


{M!-Ml)M,r^, 




)ii (^'^^^ 


^22 





(41) 



Here r^v- is the total decay rate of the right-handed neutrino A^j, 



in 



(42) 
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The CP- violating asymmetry is enhanced when the mass difference between two 
heavy right-handed neutrinos is small, although not smaller than the decay width. 

The total CP asymmetry e has an involved dependence on the complete structure of 
the neutrino matrices hy and M. However, let us assume that the lepton asymmetry is 
generated only at the decay of the lightest right-handed neutrino A^^i. This hypothesis is 
satisfied if the A^^i interactions are in equilibrium at the time of the A^2,3 decay (erasing 
any produced asymmetry), or if A'"2,3 are too heavy to be produced after inflation. As 
will be made clear in the following, this is a very plausible working assumption. In this 
case, the dynamics of leptogenesis can be described in terms of only 3 parameters : e. 
Ml, and 

m^ = [Khl)^^{Hf/M^. (43) 

The parameter mi, which determines the relevant interactions of A'^i, coincides with the 
light neutrino mass m,^^ only in the limit of small mixing angles, see eq. (PT]). 

Here we will be mostly concerned with the production mechanisms of A^i in the 
early cosmology. We will discuss a variety of these mechanisms and identify, in the 
different ranges of Mi and mi, the size of e required to generate the appropriate baryon 
asymmetry, hb/s ~ (2-9) x 10~^^. Our results can be used to check if specific particle- 
physics models for neutrino mass matrices are compatible with the various leptogenesis 
mechanisms. 

3.1 Thermal Production 

We start by considering the case in which the right-handed neutrino A^i reaches ther- 
mal equilibrium by scattering with the bath after the infiaton decay. The amount of 
lepton asymmetry generated by the N decay can then be computed by integrating the 
appropriate Boltzmann equations [^, . 

A measure of the efficiency for producing the asymmetry is given by the ratio K of 
the thermal average of the A^i decay rate and the Hubble parameter at the temperature 
T = Mi, 



2H 



(44) 



T=Mi 2xlO-3eV 

^The other neutrino mass parameters come into play only for very large values of nrii = 
(/ii//ij,)^^ {H)^/Mi {i = 2,3), when the lepton-number violating interactions mediated by N2 or A^3 
can partially erase the lepton asymmetry. 
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Here we have expressed the Ni decay width, see eq. (^), in terms of the parameters 
mi and Mi as 

Tiv, = ^^iMl (45) 

For mi ^ 2 X 10"'^ eV, K is less than unity and the decay process is out of equi- 
hbrium when Ni becomes non-relativistic. Under these conditions, the leptogenesis 
becomes very efficient. Indeed, the produced baryon asymmetry approaches its the- 
oretical maximum value obtained by assuming that each A'^i in thermal equilibrium 
eventually generates e baryons, 

nB\ 135 C(3) a . .._3 



e = 1 X IQ-^e. (46) 

S J max 

Here a is defined in eq. (|I]) and counts the number of degrees of freedom (for the 
standard model particle content, a = 28/79 and g^: = 427/4). 

For very small mi, K <^ 1 and A'^i decouples when it is still relativistic. At temper- 
atures T below Ml, the A^^i contribution to the energy density red-shifts like matter and 
therefore pNi/ptotai = ^i) / [Ag^T) . Eventually A''i matter-dominates the Universe at 
a temperature 

Tdom = ^ ^ 2 X lO-^Mi. (47) 
4 g* 

However this can happen only if A^^i does not decay beforehand. Since the decay tem- 
perature is 

T. ^ 0,,;^"^. ^ (^) {J^1^) 3 X 10» GeV, ,48) 

neutrino matter-domination occurs when 

mi < 4 X 10~^ eV. (49) 

Under these conditions, the bulk of the energy of the Universe is stored in the non- 
relativistic A^^i. At the time of decay, such energy density is converted into relativistic 
degrees of freedom whose temperature coincides with given in eq. (P3), 

pjv, = Mm^, = ^9*Tt (50) 

This yields the following baryon asymmetry 

riB riM^ 3T^ ( mi 

— = ea — - = ea—— = — — 8 x 10 e. (51) 

s s 4Mi VlO-6eV/ ^ ^ 
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In order not to reintroduce the cosmological gravitino problem, discussed in the intro- 
duction, one has to require that the temperature T^, after the right-handed neutrino 
decay is less than the maximum value allowed Tjih- This implies 

"^i< (^)'lO-^eV. (52) 

Therefore, the leptogenesis process is very efficient also when Tjih is low enough to 
suppress gravitino production. 

For K ^ 1, the departure from thermal equilibrium is reduced, and leptogenesis 
is less efficient. Larger values of e are now required. However, for mi ~ 10^^ eV, 
values of e of about 10~^ are sufficient to generate the appropriate baryon asymmetry. 
Realistic neutrino mass matrices can comfortably reproduce such values of e. Notice 
that the prediction for the baryon asymmetry depends weakly on Mi, as long as mi is 
not too large [^. This is because both the production and decay thermal rates of A''i 



at T = Ml depend only on mi, while the Mi dependence arises from lepton- violating 
H-i scattering. 

Let us now turn to discuss how primordial thermal equilibrium of A^i can be achieved. 
The first necessary condition is 

Ml < TnH. (53) 

This can be quite constraining, especially in view of the bound derived from the disrup- 
tive gravitino effects on nucleosynthesis discussed mentioned in the introduction. When 
the inequality (^) is not satisfied, one expects the number density of the A'^i-particles 
generated during the reheating stage to be quite small, making this case marginal, as 
far as the generation of baryon number is concerned. We would only like to mention 
here that such a number density depends upon the fine details of the dynamics of the 
reheating stage itself. In particular, the reheat temperature TptH is not the maximum 
temperature obtained after infiation; the maximum temperature is, in fact, much larger 
than TjiH [^. As a result, the abundance of massive particles may be suppressed only 
by powers of the mass over the temperature, and not exponentially. 

A second condition for A'^i thermalization is derived from the requirement that in- 
verse decay or production processes of the kind £g(3) Nit (mediated by Higgs-boson 
exchange) are in thermal equilibrium before A'^i becomes non-relativistic. This implies 

mi > 10-3 eV. (54) 
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This condition excludes the possibihty of the most efficient leptogenesis with K < 1. 



However, even if mi is somewhat smaller than the value indicated by eq. (|5^), a sufficient 
number of A^^i can be produced. Indeed, for mi > 10^^ eV, values of e ^ 10^^ can 
give rise to the observed baryon asymmetry. In the case of supersymmetric models, 
the constraint can be even less stringent and values e > 10~^ are sufficient for 

mi > 10-^ eV. 

The constraint on mi from thermalization can be evaded if new interactions, different 
from the ordinary Yukawa forces, bring A^i in thermal equilibrium at high temperatures. 
For instance, one could use the extra U{1) gauge interactions included in 5*0(10) GUTs. 
These interactions can produce a thermal population of A^^i if, at T = T^h, 

Tiff ^Z'^ NN) = > H. (55) 

This requires that the mass of the extra gauge boson Mz' should be close to Tjih and 
significantly lower than the GUT scale, 

*^-<(ira) 4xl«"G=V. (56) 
3.2 Production at Reheating 

Since it is very likely that the short period of preheating does not fully extract all 
of the energy density from the infiaton field, the Universe will enter a long period of 
matter domination after preheating where the dominant contribution to the energy 
density of the Universe is provided by the residual small amplitude oscillations of the 
classical infiaton field and/or by the infiaton quanta produced during the back- reaction 
processes. This period will end when the age of the Universe becomes of the order of 
the perturbative lifetime of the infiaton field. At this point the Universe will go through 
a period of reheating with a reheat temperature T^h given by the perturbative result 
in eq. (1). 

Let us suppose that the infiaton couples to A^i, either directly or through exchange 
of other particles. In this case, the infiaton decay process can generate a right-handed 
neutrino primordial population. The condition in eq. (|53|) is replaced by the weaker 
constraint 

Ml < m^, (57) 

where m^ is the infiaton mass. 
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The fate of the right-handed neutrinos produced by the inflaton decay depends on 
the parameter choice. If Mi < T^h and mi > 10~^ eV, the Yukawa couphngs are strong 
enough to bring A^^i into thermal equihbrium, and leptogenesis can proceed as in the 
usual scenario described in sect. pTT| . 

Let us now assume that the Yukawa couplings are much smaller, and that the 
right-handed neutrino decay temperature in eq. (^) satisfies < Trh, i.e. nii < 
[Trh /MiYlO'^ eV. After reheating, the A^^i behave like frozen-out, non-thermal, rel- 
ativistic particles with typical energy E^^ ~ ^^12. The A^i population will become 
non-relativistic at a temperature T^r = TrhMi/ E^^. At this moment, the energy of 
the Universe is shared between the radiation and the A^^i component, with a ratio of the 
corresponding energy densities which has remained constant between Trh and Tjsir, 



PR 



T=T„^ PR 



(5^ 

T=Tn„ 



Pr\t=Trh ^ ^9*Trh PnAt=Trh ^ ^N^^Ni = —B^n^. (59) 

Here is the inflaton number density just before decay, obtained by requiring energy 
conservation 

''^'^=^- = 30m^(l-%/2)' 
and describes the average number of Ni produced in a decay. Below T/v/j, the A^i 
density red-shifts like matter and eventually dominates the Universe at a temperature 

'''^ = (1 -%/2) (^) ^''^ 

Therefore, if 



^ ( Bt \' I Tr,H \^( 10"GeV \\ , 



then T^, > Tdom and A^i decays before dominating. In this case, the baryon asymmetry 
is determined to be 

ns 3 ea Bs Trh 

— = ea = — . 63 

s s A{l-B^/2)m^ ^ ' 

If the inequality (^) is not satisfied, A^^i matter-dominates the Universe and we recover 
the baryon asymmetry result in eq. (|5lD. 
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A necessary condition to be satisfied is that lepton-number violating interactions 
mediated by A^j (z = 1, 2, 3) exchange are out of equilibrium at the temperature of Ni 
decay, 



4 



Tal = ^G'pmiT' <H at T = T„ (64) 
where T^, is given in eq. (|4^). This implies 

/iQi^ GeV\'/'_ , /l0i2GeV\V2xl0-3eVn' . ^ 

^1 < 77 0.1 eV and mi < ^ 2 eV. 65 

Finally, we discuss the case > Tjih, i-e. mi > (Tr/^/Mi)^10~^ eV, in which A^i 
decays immediately after it is produced. In this case, the baryon asymmetry is still 
given by eq. (|63D , but the out-of-equilibrium condition of lepton- violating interactions 
has to be imposed at T = Trh- Therefore, eq. ( |65D is replaced by 

m, < ( j 3eV, z = l,2,3. (66) 

The combination of the bounds shows that lepton-violating interactions do not give 
severe constraints on the parameters. 

3.3 Production at Preheating 

As we have seen in sect. 2, right-handed neutrinos are efficiently produced in a non- 
thermal state during the preheating stage. In our numerical studies we have tacitly 
assumed that the superheavy right-handed neutrinos were stable. Of course, the para- 
metric resonance is affected by a nonvanishing decay width of the A^i. However, contrary 
to what happens for bosons where the presence of a large decay width removes the par- 



ticles from the resonance bands rendering the preheating less efficient , for fermions 
the presence of a decay width might be even beneficial. Indeed, for stable right-handed 
neutrinos the distribution function n{k) is rapidly saturated to unity and further par- 
ticle production is Pauli-blocked. However, if the decay width is large enough, the 
right-handed neutrinos may be produced at each inflaton oscillation when m{t) ~ and 
then decay right away. This will give rise to a certain amount of lepton asymmetry at 
each inflaton oscillation until the condition in eq. (^) is met; the lepton asymmetry 
would be generated in a cumulative way. Strictly speaking, however, the numerical 
calculation of fermion preheating presented in sect. 2 applies only to the case in which 
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the right-handed neutrinos have a decay hfetime larger than the typical time-scale of 
the infiaton oscillation m^^ 

The right-handed neutrinos produced during preheating may annihilate into infia- 
ton quanta. This back-reaction will render the final right-handed neutrino abundance 
smaller and therefore leptogenesis more difficult. Imposing that the back- reaction is 
negligible requires Ta ~ n^-^cTA ~ rri^, where cr^ ~ g'^/^AirMi)^. Since the energy spec- 
trum of the right-handed neutrinos is dominated by the maximum momentum generated 
at the last infiaton oscillation, we assume that the number density of right-handed neu- 
trinos is equal to the freeze-out value riiy-^ ~ Mim^/Gvr^, see eq. (P^D, and we obtain 

TA<m, ? < 10^^ (Y0^4iv) • ^^^^ 

One should also be sure that the number density of the right-handed neutrinos is 
not depleted by self-annihilations before they decay. This requires 



r. < r„. ^ , < 10- ^^^j . (69) 

Suppose now that the right-handed neutrinos A^^i decay before the infiaton energy 
density is transformed into radiation by perturbative processes. This occurs when Fat-^ > 
r^, which implies 



1/2/ Ml 



mi > i^^j 1 X 10-^ eV. (70) 

A crucial point is that, after the generation of non-thermal right-handed neutrinos 
at the preheating stage, the ratio of the energy densities of A^^i and infiaton quanta 
remains constant: Pni/p^} = {pni/ P<i>)ph, where {pni/ P<f>)ph is the ratio generated at the 
preheating stage. Since the energy density of the Universe is dominated by infiaton 
oscillations, after preheating we obtain 

V P0 Jph 8^ 

Our numerical results discussed in sect. 2 indicate that, for very large q and Mi, 

(Pjh\ ^lo-iig for Mi<^m<^. (72) 
XP^Jph ^ 
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At tTVi ~ ^nI right-handed neutrinos decay and the energy density piVi is converted 
into a thermal bath with temperature 



—g^T = pN,, (73) 
where is computed at -ff = Ftvi- Using eqs. ([7l| ) and (|72j), we find 

Before the infiaton decay, this thermal bath never dominates the energy density of the 
Universe since p^-^ <^ at ~ F ^r-^ and the energy density in the infiaton field is red- 
shifted away more slowly than the radiation. Notice that at this time the asymmetry 
is still in the form of lepton number since the sphalerons which are responsible for 
converting the lepton asymmetry into baryon asymmetry are still out-of-equilibrium 
at T = T. One might be worried that, since T is usually larger than Trh? too many 
gravitinos are produced at the stage of thermalization of the decay products of the 
right-handed neutrino. However, one can estimate the ratio 113/2/3 after reheating to 
be of the order of 10~^^(g/10^°)^/^ (Tr///10^° GeV), which is quite safe. However, we 
have to require that the lepton-number violating processes within the thermal bath at 
temperature T are out-of-equilibrium in order not to wash out the lepton asymmetry 
generated by the right-handed neutrino decays. Therefore, we demand that 

^Glm^T^ <H at T = T. (75) 



AL 



Tf' 



This implies 



/lOio\3/2. 15 y.^xlO-^eVn' , 
mi < — ^ 6 X 10"^ eV. 77 

At H < FjVi, the ratio riL/n^ keeps constant until the time ~ F^""^ when the 
infiaton decays and the energy density in the infiaton field P4i{t,i)) is transferred to 
radiation. After reheating we obtain the following lepton asymmetry to entropy density 
ratio 

nL_ _ 'iriL Trh _ 3 e Trh ( Pni \ 



s 



Ap^{t^) A Ml \p^)^^ 
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Using eq. (|72D, the baryon asymmetry can be expressed as 
riB ( Trh 



s VlQiOGeV 



\ /lO^s GeV\ f q \ 7 



(79) 



The result in eq. (|79]) shows that the preheating production mechanism can lead to 
a successful leptogenesis even for Mi as large as 10^^ GeV, if eq ~ 10^. Values of q as 
large as 10^'' correspond to a perturbative coupling between A^^i and the infiaton g"^ ~ 0.4 
and are compatible with the constraints in eqs. (|68D and (|69l). Values of e of the order 
of 10~^ are quite large, but can be attained with realistic neutrino mass matrices. For 
values of Mi so close to the GUT scale, we expect that all the right-handed neutrino 
Majorana masses are comparable in size. Moreover, the atmospheric neutrino results, 
together with the requirement of perturbative Yukawa couplings, indicate that at least 
one Majorana mass is less than about 8 x 10^^ GeV. This situation of comparable 
Majorana masses and some large Yukawa couplings naturally leads to large values of 
e. Also, notice that the out-of-equilibrium condition is automatically satisfied in the 
preheating scenario for all 3 right-handed neutrinos, in a large range of parameters. 

Let us finally consider the case in which A'^i decays after the reheating process and 
the inequality ([70| ) is not satisfied. Again, it is important to establish whether A'^i 
dominates the Universe before decaying. Since the infiaton energy density is converted 
into radiation and the A'^i are non-relativistic, the temperature at which p^^ dominates 
is given by 

Tdom = ( ^— ) Trh- (80) 
V P't' Jph 

Therefore, for > Tdom, for 

m.>(^y(^yiO''eV, (81) 



Ml ; vioio. 

the estimate for the baryon asymmetry in eq. (|79|) is still valid. Otherwise, for < Tdom, 
we obtain the result in eq. (151]) . Notice that, in this case, the constraint T^ < Trh is 
automatically satisfied. 

3.4 Comparison of the Different Production Mechanisms 

We want to compare here the different mechanisms discussed in this section for lep- 
togenesis from A'^i decay. We summarize their most important features and estimate 
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the size of the CP-violating parameter e necessary to reproduce the observed baryon 
asymmetry. 

Thermal production. This is the right-handed neutrino production mechanism 
usually considered in the literature for conventional leptogenesis. In the range 10~^ eV < 
mi < 10~^ eV (or 10~^ eV < mi < 10~^ eV in the case of the minimal supersymmetric 
model) and for Mi < T^h, the thermal production of unstable A^^i is efficient, and 
values of e in the range lO^'^-lO^^ can account for the present baryon asymmetry. The 
boundaries of the allowed region of neutrino mass parameters are determined as follows. 
For small values of mi, the A^^i production rate is suppressed and larger values of e are 
required. For large mi, the Yukawa couplings maintain the relevant processes in thermal 
equilibrium for longer times, partially erasing the produced asymmetry. The values of 
Ml are limited by the reheat temperature after inflation Tj^h, which in turn is bounded 
by cosmological gravitino considerations to be below 10^-10^° GeV. 

Production at reheating. If A'^i is directly or indirectly coupled to the infia- 
ton, the decay of the small amplitude oscillations of the classical infiaton field at the 
time of reheating can produce a right-handed neutrino population. This production 
mechanism enables us to extend the leptogenesis-allowed region to neutrino mass pa- 
rameters which correspond to non-thermal A^i populations. In particular. Mi can 
be as large as m^ 10^^ GeV. The observed baryon asymmetry is reproduced for 
e ~ 10^^(10^° GeV/TR//)(lO^^/i?0), where is the average number of A''i produced 
by a single infiaton decay. 

Production at preheating. The non-perturbative decay of large infiaton oscilla- 
tions during the preheating stage can produce a non-thermal population of very massive 
right-handed neutrinos. In this case, the range of Mi can be extended to values close 
to the GUT scale, while mi is bounded by the condition that lepton-number violating 
interactions are out of equilibrium after the A''i decay, see eq. ([77|). A successful lepto- 
genesis requires e ~ 10-^(10^° GeV /Trh) {Mi/ 10^^ GeV){10^° /q), where q is related to 
the initial infiaton configuration and is defined in eq. (]23|). 

In conclusion, leptogenesis provides an interesting and simple way to explain the 
present cosmic baryon asymmetry. The study of the different mechanisms in which it 
can be realized provides us with precious information on the neutrino mass parameters 
and the early history of the Universe. 
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